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ABSTRACT 

The size-luminosity relation of early-type Brightest Cluster Galaxies (BCGs), R e oc L - 88 , is steeper 
than that for the bulk of the early-type galaxy population, for which R e oc L 68 . This is true if quantities 
derived from either deVaucouleur or Sersic fits to the surface brightness profiles are used. Contamination 
from an intra-cluster light component centered on the BCG, with similar parameters to what has been 
seen in some recent studies, is not able to account for this difference. In addition, although BCGs are 
hardly offset from the Fundamental Plane defined by the bulk of the early-type population, they show 
considerably smaller scatter. The larger than expected sizes of BCGs, and the increased homogeneity, 
are qualitatively consistent with models which seek to explain the colors of the most massive galaxies 
by invoking dry dissipationless mergers, since dissipation tends to reduce the sizes of galaxies, and wet 
mergers which result in star formation would tend to increase the scatter in luminosity at fixed size and 
velocity dispersion. Furthermore, BCGs define the same g — r color-magnitude relation as the bulk of the 
early-type population. If BCGs formed from dry mergers, then BCG progenitors must have been red for 
their magnitudes, suggesting that they hosted older stellar populations than typical for their luminosities. 
Our findings have two other consequences. First, the R e ~ L relation of the early-type galaxy population 
as a whole (i.e., normal plus BCG) exhibits some curvature: the most luminous galaxies tend to have 
larger sizes than expected from the R e oc L om scaling — some of this curvature must be a consequence 
of the fact that an increasing fraction of the most luminous galaxies are BCGs. The second consequence 
is suggested by the fact that, despite following a steeper size- luminosity relation, BCGs tend to define a 
tight relation between dynamical mass R e a 2 /G and luminosity. Although this relation is slightly different 
than that defined by the bulk of the population, the fact that their sizes are large for their luminosities 
suggests that their velocity dispersions are small. We find that, indeed, BCGs define a shallower a — L 
relation than the bulk of the early-type galaxy population. This shallower relation suggests there may 
be curvature in the correlation between black hole mass and velocity dispersion; simple extrapolation of 
a single power law M. — a relation to large a will underestimate M m . 

Subject headings: galaxies: elliptical — galaxies: evolution — galaxies: fundamental parameters — 
galaxies: photometry — galaxies: stellar content 



1. INTRODUCTION 

Understanding why massive early-type galaxies are red 
and dead has proved to be difficult; the problem is to ar- 
range for star formation to occur at higher redshift than 
the actual assembly of the stars into a single massive 
galaxy. The most recent galaxy formation models arrange 
for this to happen by a combination of two processes: dry 
mergers and AGN feedback (Hopkins et al. 2005; Croton 
et al. 2005; De Lucia et al. 2005; Robinson et al. 2006; 
Bower et al. 2006). The dry merger hypothesis assumes 
that the assembly of massive galaxies occurs by merging 
smaller progenitor systems of old stars without additional 
star formation (which would otherwise lead to bluer col- 
ors). This happens either because the merging units were 
themselves gas poor, or because AGN feedback has pre- 
vented the hot gas reservoirs of the progenitors from cool- 
ing and forming stars after the merger. Together, these 
processes allow massive galaxies to be built from smaller 
systems while still remaining red. 
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The dry merger hypothesis is most necessary for the 
most massive galaxies. These tend to be the brightest 
galaxies in clusters (BCGs). Our primary focus will be 
the sizes of these BCGs, although we also study a num- 
ber of other scaling relations. This is because, in a 'wet' 
merger, gas from the merging components is able to dissi- 
pate energy, cool, contract and eventually make new stars. 
Dry mergers on the other hand have no energy loss mech- 
anism (other than dynamical friction) ; having suffered less 
dissipation, their stellar components are not as centrally 
contracted, so the optical sizes of dry merger products 
are expected to be larger than if the mergers were wet 
(e.g. Kormendy 1989; Capelato et al. 1995; Nipoti et 
al. 2003), although the sizes may also reflect the orbital 
parameters of the mergers which formed the BCG (e.g. 
Boylan-Kolchin et al. 2006). 

The structural properties of BCGs have been the object 
of much previous work. Compared to normal early type 
galaxies, BCGs have larger than expected radii for their 
surface brightnesses (Thuan & Romanishin 1981; Hoes- 
sel et al. 1987; Schombert 1987, 1988; Oegerle & Hoes- 
sel 1991), and smaller than expected velocity dispersions 
for their luminosities (Malumuth & Kirshner 1981, 1985; 
Oegerle & Hoessel 1991). Crawford et al. (1999) describe 
the results of an extensive study of the spectra of X-ray 
selected BCGs, and Brough et al. (2005) find that BCGs 
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Fig. 1. — Images (top) and surface brightness profiles (bottom) for a set of C4 BCGs which are well fit by a single 
deVaucouleur profile. The fit is good in both the g and r bands. Horizontal dashed lines show one percent of the sky 
brightness in the two bands. Bottom half of bottom panels shows the associated g — r color gradient — the gradient is 
expected to be weak if the formation was dominated by mergers. 



of X-ray luminous clusters tend to be have larger than ex- 
pected sizes. With the exception of Crawford et al., these 
studies were based on sample sizes of order 50 or smaller, 
so it is interesting to ask if the larger samples now available 
show the same trends. 

We will study the objects which are identified as BCGs 
in the C4 cluster catalog of Miller et al. (2005). This cata- 
log is complete out to z — 0.12, contains about 700 BCGs, 
and will contain about 1500 BCGs by Fall 2006. Section 2 
describes the sample and provides a brief summary of our 
fits to the photometric properties of BCGs and normal 
early-type galaxies. Details are provided in Hyde et al. 
(2006). These new fits were necessary, as the SDSS photo- 
metric reductions of objects in crowded fields suffer from 
sky subtraction problems (e.g. Hyde et al. 2006, Lauer et 
al. 2006). In Section 3 we compare a variety of BCG scal- 
ing relations to those of the bulk of the early-type galaxy 
population. This section shows that BCGs define a steeper 
size-luminosity relation than does the bulk of the early- 
type galaxy population; that although BCGs have larger 
than expected sizes, they still define a tight relation be- 
tween luminosity and dynamical mass Ra 2 /G, and so the 
a — L relation of BCGs curves away from the power-law 
scaling defined by the bulk of the early-type galaxy pop- 
ulation. It also shows the location of BCGs in the same 
Fundamental Plane defined by normal early-type galaxies, 
and shows that they define the same color-magnitude re- 
lation as normal early- type galaxies. Appendix A presents 
a comparison with the color-magnitude relation predicted 
by the semi-analytic galaxy formation models of Croton et 



al. (2006) and Bower et al. (2006). 

Section 4 studies two systematic errors which might have 
caused the steeper R e — L relation we see, and concludes 
that they are probably not to blame. The first possibil- 
ity is that the deVaucouleur fit is systematically worse for 
BCGs than it is for the bulk of the population. The sec- 
ond possibility is more subtle. Galaxy clusters contain a 
substantial population of stars which are not associated 
with a galaxy (e.g. Gonzalez et al. 2005; Zibetti et al. 
2005). So it is not unreasonable to ask if this intracluster 
light is contaminating our estimates of the sizes and lumi- 
nosities of BCGs, and hence the inferred R e — L and a — L 
relations. 

At the very largest luminosities (M r < —23.5), the R e — 
L relation of the early-type galaxy population (normal + 
BCG) appears to curve upwards toward larger sizes from 
the R e oc L ' 68 scaling. If an increasing fraction of these 
objects are BCGs, then this curvature is not unexpected 
in view of our finding that BCGs follow a steeper size- 
luminosity relation. Appendix B presents evidence that, 
at large luminosities, an increasing fraction of galaxies are 
indeed BCGs. However, our results do not exclude the 
possibility that there is, in addition, intrinsic curvature in 
the R e — L relation of early-types which are not BCGs. 

A final section summarizes our findings and discusses 
the implications for galaxy formation models, and for 
studies which seek to estimate the mass of a supermas- 
sive black hole from the stellar velocity dispersion of the 
bulge which surrounds it. Where necessary, we assume a 
flat cosmological model with SIq = 0-3, A = 1 — fio and 
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Fig. 2. — Images (top) and surface brightness profiles (bottom) for a set of C4 BCGs which are not well fit by a 
deVaucouleur profile. The deviation is strongest at large radii (faint surface brightness). Although Sersic profiles provide 
better fits, Gonzalez et al. (2006) show that the sum of two deVaucouleur components generally provide an even better 
description of most such profiles. 



H = 70h 70 km s^Mpc" 1 . 
2. SAMPLE SELECTION AND FITS TO LIGHT PROFILES 

We have selected the objects identified as BCGs in the 
C4 cluster catalog of Miller et al. (2005). This catalog 
was based on the Second Data Release of the SDSS: DR2. 
Of the 748 BCGs in the catalog, spectra are available for 
580 (fiber collisions mean some BCGs were missed). Some 
of these are duplicates (i.e. some clusters were listed more 
than once), and the number which are matched to a par- 
ent catalog based on SDSS DR4 is 403. Of these, a vi- 
sual inspection of the fields centred on C4 BCGs fainter 
than M r < —21 often showed a brighter galaxy nearby, 
especially when the C4 BCG was a spiral. If we restrict 
attention to objects with M r < —21, then 286 have SDSS 
reported velocity dispersions, and have at least 10 other lu- 
minous galaxies (M r < -21) within l/i^Mpc. The SDSS 
does not estimate velocity dispersions from spectra which 
have strong emission lines, so our requirement that the 
SDSS reports a velocity dispersion means that the objects 
we select are almost certainly not AGN. Removing the fi- 
nal cut on 'richness' increases the sample of BCGs to 324, 
but does not change the results which follow. 

Visual inspection of the SDSS images of these objects 
shows that 25% actually have spiral arms; the absence of 
emission lines in the spectra of these objects is presumably 
because the SDSS fiber, with a radius of 1.5 arcsec, only 
sees the inner bulge. In what follows, we will be careful 
to distinguish the scaling relations of the ~ 215 genuine 
early- type BCGs from the rest. 



It happens that BCGs are a class of object for which 
the SDSS photometric reductions are particularly unreli- 
able: the reported magnitudes are in error by far more 
than the quoted 0.02 mags, with some discrepancies as 
large as 1 mag (Hyde et al. 2006). This is due to a com- 
bination of factors: because BCGs are often in crowded 
fields, the SDSS reductions tend to overestimate the sky 
level. This results in underestimates of the BCG magni- 
tude and half-light radius, which are particularly severe for 
BCGs with large half-light radii. In addition, the surface 
brightness distributions of BCGs do not always follow pure 
single component deVaucouleurs profiles, which the SDSS 
reductions assume. For these reasons, we have performed 
our own photometric reductions, and examined the image 
of each BCG by eye. A more detailed discussion of the 
various issues involved is provided by Hyde et al. (2006), 
whose results we use below. 

Because we will be interested in comparing the proper- 
ties of BCGs to those of other early-type galaxies, we have 
also reanalyzed the photometric properties of the sample 
of 9000 SDSS early-type galaxies defined in Bernardi et 
al. (2003a). It is these parameters (i.e. those based on 
Hyde et al. 2006), rather than those output from the 
SDSS pipeline, which we use in what follows. On average, 
the Hyde et al. magnitudes are 0.1 mags brighter than 
those reported by the SDSS DR4 pipeline, and the sizes 
are about 10% smaller. In crowded fields, and especially 
for bright objects, Hyde et al. report both brighter magni- 
tudes and larger sizes than does the SDSS. When compar- 
ing the scaling relations of these objects with those for the 
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C4 BCGs, we restrict the redshift range to that over which 
the C4 catalog is complete (z < 0.12; except for the Fun- 
damental Plane, this restriction makes little difference). 

The most important finding of this reanalysis of BCG 
photometry is that there appear to be two classes of early- 
type BCG light profiles. Examples of these two types are 
shown in Figures 1 and 2. The first type is well fit by a 
single component deVaucouleurs profile, whereas the sec- 
ond is not. This second type of profile is similar to that 
associated with cD galaxies. Section 4 discusses these two 
profile types in more detail. Meantime, we will be careful 
to treat these two populations separately in what follows. 

3. BCG SCALING RELATIONS 
3.1. The size-luminosity relation of BCGs 

Figure 3 shows the correlation between luminosity and 
half-light radius for the C4 BCGs: diamonds represent 
objects with spiral morphologies, and filled circles repre- 
sent E/SOs. Solid line shows the R e — L relation defined 
by the bulk of the early-type galaxy population — BCGs 
clearly define a steeper relation. (The steeper relation is 
also seen if we use the original SDSS photometric reduc- 
tions which we believe to be incorrect — see bottom right 
panel of Figure B13.) In what follows, we focus on the 
early-type BCGs, but note that the spiral BCGs define an 
B e — L relation that is at least as steep. 

Figure 4 shows the R e — L relation for the early-type 
BCGs, now subdivided by whether (filled circles) or not 
(open triangles) the surface brightness profile is well fit by 
a single deVaucouleurs profile. Jagged solid line shows the 
median BCG size (whether or not it is well-fit by a deVau- 
couleur profile) in a few bins in luminosity, and dashed 




-20 -21 -22 -23 -24 -25 
M r [mag] 



Fig. 3. — C4 BCGs define a steeper size-luminosity rela- 
tion (symbols) than does the bulk of the early-type galaxy 
population (solid line). Open diamonds and filled circles 
represent BCGs with spiral and E/SO morphologies, re- 
spectively. In what follows, we only consider the E/SO 
objects. 



line shows a power-law fit to the size-luminosity relation 
for these BCGs: 

^log 10 .R e |M r \ = -0.885 (M r + 21)/2.5 + 0.230. (1) 

This fit should be compared to the straight solid line which 
shows the scaling that provides a good description of the 
bulk of the early-type galaxy population: 

(log 10 R e \M r \ = -0.681 (M r + 21)/2.5 + 0.343. (2) 

Clearly, BCGs define a steeper relation — they have larger 
sizes than do normal early-type galaxies of similar lumi- 
nosity. This trend is consistent with the hypothesis that 
BCGs have formation histories with less dissipation than 
the bulk of the early-type population. 

The fits above are based on the photometric reductions 
described by Hyde et al. (2006). Using the SDSS photom- 
etry (deVaucouleur magnitudes and sizes) instead gives 

log 10 i? e |M r ^) = -0.764 (M r + 21)/2.5 + 0.395 

for the BCGs and 

Aog 10 R e \M r ^ = -0.623 {M r + 21)/2.5 - 0.402 

for the total early-type galaxy population. 
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Fig. 4. — Early-type C4 BCGs define a steeper size- 
luminosity relation (symbols) than does the bulk of the 
early- type galaxy population (solid line). Filled circles rep- 
resent objects which are well-described by a deVaucouleurs 
profile, and open triangles represent objects which are not. 
Dashed line shows a fit to the R e —L relation using all these 
BCGs (i.e. both filled and open symbols), and jagged solid 
line shows the median BCG size as a function of luminos- 
ity. 
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3.2. The mass-luminosity relation 

The larger than expected sizes of BCGs have interesting 
consequences. If these objects are virialized, then one ex- 
pects a relatively tight correlation between luminosity and 
the dynamical mass cx R e a 2 /G. Figure 5 shows that this 
correlation is indeed rather tight — the inset shows that 
BCGs define a tighter relation than does the bulk of the 
early-type galaxy population. This smaller scatter is also 
consistent with dry-merger formation histories; one might 
have expected 'wet' mergers to result in star formation, 
which would increase the scatter in luminosities at fixed 
dynamical mass. 

Figure 5 also shows that BCGs define a slightly different 
correlation between luminosity and R e a 2 than does the 
total early-type galaxy population. The actual relations 
are 



log 1Q R e a 2 \M r ) = -1.114 (M r + 21)/2.5 + 4.822 (3) 



with a scatter around the relation of 0.11 dex for the BCGs 
and 



log 10 R e a z \M 7 



-1.185 (M r + 21)/2.5- 4.697 (4) 



with a scatter around the relation of 0.34 dex for the total 
early-type galaxy population. 

3.3. The velocity dispersion-luminosity relation 
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Fig. 6.— The a — L relation of C4 BCGs (filled and 
open symbols represent objects which are and are not well 
described by a deVaucouleurs profile) appears to flatten 
at high luminosities. This flattening may be related to 
studies of curvature in the M, — a relation. Solid and 
dashed lines show linear fits to this relation for the early- 
type galaxy population and for this BCG sample. Solid 
jagged lines show the median value as a function of L for 
objects which are (lower line) and are not well-fit (upper 
line) by a single deVaucouleur profile. 



Fig. 5. — Correlation between dynamical mass R e o~ 2 and 
luminosity in the bulk of the early-type population (solid 
line) and in the BCG sample (dashed line); filled circles 
and open triangles represent BCGs which are and are not 
well-fit by a deVaucouleurs profile. Solid jagged line shows 
the median value as a function of L. BCGs define a slightly 
shifted and even tighter correlation between dynamical 
mass R e cr 2 and luminosity than does the bulk of the early 
type galaxy population (inset). 
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If the R e — L relation for BCGs is steeper than for the 
bulk of the early-type galaxy population, then the tight 
correlation between L and dynamical mass shown in Fig- 
ure 5 suggests that the a— L relation for BCGs will be shal- 
lower. Figure 6 shows that this does appear to be the case. 
Once again, the solid line shows the relation for the full 
sample, and the dashed line shows a fit to the BCGs (filled 
and open symbols). A running average over small bins 
in luminosity shows evidence for curvature in the BCG 
a — L relation, in the sense that velocity dispersion in- 
creases only weakly with luminosity once a >250 km/s. 
However, notice that objects which are not well- fit by a 
single deVaucouleur profile appear to have slightly larger 
velocity dispersions. 

3.4. The Fundamental Plane 

Numerical simulations of a dry merger scenario sug- 
gest that the final object will remain on the Fundamental 
Plane defined by the bulk of the early-type population (e.g. 
Boylan-Kolchin et al. 2006). Figure 7 shows the location 
of the C4 BCGs in the Fundamental Plane defined by the 
bulk of the early-type galaxy population: 

log 10 R e = 1.307 a + 0.763 p, r /2.5 - 8.345. (5) 

This is the "orthogonal" Fundamental Plane and it differs 
from the result presented in Bernardi et al. (2003b) mainly 
because we have restricted our sample to z < 0.12 (the 
new photometry changed the zero-point but not the slope 
of the relation) . 

The BCGs which are well fit by a single deVaucouleur 
profile lie on this same Fundamental Plane, whereas the 
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Fig. 7. — Location of C4 BCGs (filled and open symbols 
as in previous figures) with respect to the Fundamental 
Plane defined by the bulk of the early-type galaxy popula- 
tion (dots). Inset shows that BCGs are distributed more 
tightly around this FP than the bulk of the early-type 
galaxy population: the rms scatter is 0.04 dex compared 
to 0.058 dex. 



others are slightly offset (—0.025 dex) from it. The inset 
shows that the scatter around this Plane is significantly 
reduced — the scatter orthogonal to the Plane decreases 
from 0.058 dex (solid line), to 0.046 dex (dotted line), to 
0.04 dex (dashed-line) if only BCGs which are well fit by 
a deVaucouleur profile are used. This reduced scatter for 
BCGs is consistent with a dry merger formation history — 
even a small amount of star formation associated with a 
wet merger has the potential to change the surface bright- 
nesses significantly, resulting in increased scatter about the 
Plane. 

3.5. The color-magnitude relation 

Early-type galaxies define a tight correlation between 
color and luminosity. Figure 8 shows that BCGs follow 
the same correlation, though with slightly smaller scat- 
ter. The Appendix presents a comparison with the color- 
magnitude relation predicted by semi-analytic galaxy for- 
mation models. The models suggest that BCGs should lie 
slightly blueward of the relation defined by the bulk of the 
population — an effect which is not seen in Figure 8. 

Whereas the luminosity of an object formed from a dry 
merger is expected to be the sum of the luminosities of 
the objects involved in the merger, the merger is not ex- 
pected to change the color (there being no associated star 
formation). Therefore, if BCGs formed from dry mergers, 
then the fact that they define the same color-magnitude 
relation as the bulk of the population suggests that the 
progenitors of BCGs were red for their luminosities. Ob- 
jects which lie redward of the color-magnitude relation are 
expected to be older (Kodama et al. 1998; Bernardi et al. 
2005). Thus, Figure 8 is not inconsistent with the expec- 
tation that the oldest stars in a halo find their way to the 
bottom of the potential well, the location of which is often 
traced by the BCG. 

Figure 8 shows that BCGs also lie on approximately 
the same er-color correlation as the bulk of the early-type 
population. This is consistent with a model in which the 
color is not altered by the merger, and the change to er 
is also weak (recall Figure 6 suggests that dry mergers do 
not change a as much as they change L). 

The color-magnitude relation for the bulk of the early- 
type galaxy population is entirely due to the color-tr and 
cr-magnitude correlations (Bernardi et al. 2005). The fig- 
ures above suggest that this is not true for BCGs. Specif- 
ically the flatter a — L relation for BCGs, combined with 
the fact that the color-magnitude and color— a relations 
show little difference from the main sample implies that 
the correlation between color and magnitude for BCGs is 
substantially stronger than expected from the BCG color- 
a and a- magnitude correlations. 

4. TESTS OF SYSTEMATIC EFFECTS 

This section studies two possible sources of systematic 
trends and concludes that they do not cause the steeper 
R e — L relation we see in our early-type BCG sample. 

4.1. Sersic profile fits 

We have considered the possibility that the steeper 
R e — L relation we see arises because the deVaucouleur 
profile provides a systematically worse description of early- 
type surface brightness profiles at large L. Figures 9-10 
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M r [mag] Log 10 a [kms '] 

Fig. 8. — Location of C4 BCGs (filled and open symbols represent objects which are and are not well-fit by a deVaucouleurs 
profile) with respect to the color-magnitude (left) and color-cr (right) relations defined by the bulk of the early-type galaxy 
population (solid line). Dashed line shows a linear fit of color as a function of absolute magnitude (or velocity dispersion), 
and jagged line with error bars shows the median color in a few bins in luminosity (or velocity dispersion). Inset shows 
that BCGs are not offset from the relation defined by the bulk of the population, though the relation they define is slightly 
tighter. 




Fig. 9. — Left: Correlation between Sersic parameter n and the total luminosity obtained from fitting the Sersic form to 
the surface brightness profiles. Filled and open symbols are the same as in previous figures. The filled symbols tend to 
have n = 4, the value which represents a deVaucouleur profile, whereas objects which show evidence for two components 
tend to have large n. When the deVaucouleur profile is a poor fit (open symbols), the Sersic fits result in significantly 
larger total luminosities (middle panel) and sizes (right). 
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show results obtained from fitting Sersic profiles to the full 
early-type sample, as well as to the BCGs. Figure 9 shows 
that objects which we identified as potentially having two 
components (like those shown in Figure 2) tend to have 
large values of n. Figure 9 shows that, for those objects 
which we flagged as having two components, the luminosi- 
ties and half-light radii derived from the Sersic fits tend to 
be substantially larger than when the fit is restricted to a 
deVaucouleur profile. In some cases the difference in total 
luminosity is as large as one magnitude — this systematic 
difference dwarfs the 0.02 mag uncertainty usually quoted 
for SDSS photometry. Figure 10 shows the analog of Fig- 
ure 4: the R e — L relation obtained from the Sersic profile 
fits. This illustrates that BCGs define a steeper relation 
than the bulk of the population even when parameters de- 
rived from Sersic fits are used. 

4.2. The effect of intra-cluster light 

The previous section showed that BCGs define a steeper 
size-luminosity relation than does the bulk of the early- 
type galaxy population. However, it is known that there 
is a component of the total luminosity of a cluster which is 
not associated with any particular galaxy in it: this intra- 
cluster light (ICL) tends to have a smooth profile centered 
on the cluster center (e.g. Malumuth & Kirshner 1981; Zi- 
betti et al. 2005; Gonzalez et al. 2005; Krick et al. 2006). 




Fig. 10. — The R e — L relation obtained from fitting Ser- 
sic profiles to the surface-brightness distribution (compare 
Figure 4) . Solid line shows the R e — L relation defined by 
the bulk of the early-type galaxy population and dashed 
line shows the steeper relation defined by the C4 BCGs 
when the fit is restricted to n — 4, a deVaucouleur profile; 
symbols show that C4 BCGs define a steeper relation even 
when Sersic parameters are used. Filled symbols represent 
objects which are well fit by a single component deVau- 
couleur profile in the previous figures, and open symbols 
show the other early-type BCGs. 



So it is natural to ask if this ICL profile is affecting our 
estimates of the sizes and luminosities of BCGs, thus al- 
tering the correlation between size and luminosity. This 
effect would be particularly severe for BCGs which lie at 
the cluster center. For such objects, it may be more rea- 
sonable to treat the surface brightness profile as the sum of 
two components, an inner one which represents the BCG 
itself, and an outer one which represents the ICL. 

To estimate the effect of the ICL on our results, we 
have used two results from the recent literature. Zibetti et 
al. (2005) used a stacking technique to detect the ICL 
with high significance in SDSS clusters, and concluded 
that the ICL component has rather low surface bright- 
ness (fi r > 25 mag arcsec -2 ): single SDSS exposures in 
the r— band are not sufficiently deep to detect it unam- 
biguously, which is why a stacking analysis was required. 
This suggests that contamination from the ICL is unlikely 
to be a strong effect in our case. However, recent work 
by Gonzalez et al. (2005) shows that BCG light profiles 
deviate from a deVaucouleur fit at about the same surface 
brightness values (once the difference between their / band 
and our r band has been accounted for) at which we see 
deviations (c.f. Figure 2). 

To check if the ICL is affecting our results, we made 
mock observations of BCGs of various luminosities as fol- 
lows. For a given BCG luminosity, we require the half-light 
radius to satisfy R e oc L 68 . We use this size and luminos- 
ity to generate a deVaucouleurs surface brightness profile. 
To this profile we add an ICL component, the parame- 
ters of which are motivated by the study of Zibetti et al.. 
Namely, the ICL component is modeled as a second de- 
Vaucouleurs profile, with the same integrated luminosity 
as the first component, but with half-light radius twenty 
times larger. We then fit the summed image with a single 
component deVaucouleurs profile (since this is, in effect, 
what the SDSS pipeline does), and record the best-fit val- 
ues of luminosity and half-light radius. Upon repeating 
this procedure for a range of luminosities, we can compare 
the i?BCG+iCL — ^bcg+icl relation with the input one. 

In general, the ICL component tends to increase both 
the size and the total luminosity of the BCG. Our exper- 
iments indicate that the net effect of the ICL is to move 
the BCG along the R e — L relation. Therefore, we find 
that the ICL component does not lead to a steepening of 
the R e — L relation. This conclusion remains unchanged 
if we use Sersic profiles instead, or the particular choice of 
Petrosian quantities used by the SDSS. (Because the Pet- 
rosian quantities do not make strong assumptions about 
the profile shape, one might have thought that they were 
particularly well suited to this comparison. However, the 
Petrosian quantities output by the SDSS pipeline are nei- 
ther seeing-corrected nor have the correct sky subtraction, 
making them less attractive.) 

This conclusion also remains unchanged if we assume 
that the ICL component follows a NFW profile (i.e., the 
profile traced by the dark matter in simulations of dissi- 
pationless clustering). For NFW concentrations between 
5 and 10 (the range associated with cluster mass halos), 
and £icl/£bcg between 1 and 10, the net effect of the 
ICL is to move the BCG approximately along the R e — L 
relation. At larger Xicl/^bcGj the ICL component re- 
sults in a flattening rather than a steepening of the R e — L 
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relation. 

5. DISCUSSION 

Whether or not BCGs comprise a special population 
has been the subject of considerable debate for almost 
fifty years (Scott 1957; Sandage 1976; Schechter & Pee- 
bles 1976; Tremaine & Richstone 1977; Bhavsar & Barrow 
1985; Oegerle & Hoessel 1991; Postman & Lauer 1995; 
Loh & Strauss 2006), with different authors coming to dif- 
ferent conclusions. Semi-analytic galaxy formation models 
have long made the assumption that the central object in 
a halo is indeed special (Croton ct al. 2006; Bower et al. 
2006). And recently halo model (see Cooray & Sheth 2002 
for a review) interpretations of galaxy clustering strongly 
suggest that this is indeed the case (Berlind et al. 2005; 
Skibba et al. 2006). 

One of the reasons the subject has received so much 
attention (and controversy!) is that BCGs appear to be 
quite homogeneous — being big and bright, they may be 
used as standard candles out to large distances. Their 
attractiveness as standard candles increases if their nature 
is understood, since the next generation of photometric 
galaxy surveys will identify large numbers of clusters to 
redshifts of order 0.7. This means that large BCG samples 
will also become available. If their nature is understood, 
then the BCGs themselves may provide complementary 
information about cosmological parameters. 

BCG samples are only just becoming sufficiently large 
that unambiguous conclusions about their nature can be 
drawn. The particular question we addressed is whether or 
not BCG formation is dominated by dissipation or by dry 
mergers. The dry merger hypothesis seeks to explain why 
massive early-type galaxies, which are expected to have as- 
sembled their stellar mass recently, host the oldest stellar 
populations — even when the age estimate is luminosity- 
weighted. Recent observations suggest that dry mergers 
are a common formation mechanism for massive early-type 
galaxies (e.g., van Dokkum 2005). Since most BCGs are 
early-types, it is reasonable to ask if they show some sig- 
nature of unusual formation histories. 

Our analysis suggests that the answer to this question is 
'yes'. The optical surface brightness profiles of C4 BCGs 
appear to be of two types — one has the classical deVau- 
couleur profile (Figure 1), and the other shows evidence for 
a second, more extended, lower surface brightness compo- 
nent (Figure 2). The BCGs which are well-fit by classic 
deVaucouleur profiles have larger half-light radii than nor- 
mal early- type galaxies of similar luminosities (Figure 4). 
The other BCGs also appear to follow a steep R e — L 
relation. The low surface brightness component in these 
objects may be due to light from the ICL. If so, then the 
steeper R e — L relation is also true of the BCG compo- 
nents of these objects, since the ICL tends to move ob- 
jects along the R e — L relation, rather than change its 
slope (Section 4.2). 

We also examined the R e — L relation of the normal 
early-type galaxy population, with a view to seeing if the 
most luminous galaxies tend to have abnormally large sizes 
whether or not they are BCGs. Although there is clear 
evidence for curvature in the relation defined by the full 
early-type sample (BCG + non-BCG), the evidence for 
curvature in the early-type sample once BCGs had been 



removed was much less compelling (Figure B13). 

The R e — L relation of the bulk of the early-type galaxy 
population is steeper than the correlation between half 
mass radius and mass of dark matter halos. This is consis- 
tent with models in which galaxies in low mass halos have 
suffered more dissipation (e.g. Kormendy 1989; Dekel & 
Cox 2006). If massive galaxies have suffered less dissipa- 
tion than less massive galaxies, then our results indicate 
that BCGs have suffered even less dissipation than other 
massive galaxies of similar luminosities. Although some 
of the steepening of the R e — L relation may be due to 
anisotropic mergers (e.g., Boylan-Kolchin ct al. 2006), the 
larger than expected sizes of BCGs also suggest less dissi- 
pation, so they are consistent with the dry merger forma- 
tion hypothesis. In addition, BCGs show a substantially 
smaller scatter around the Fundamental Plane defined by 
the bulk of the early-type population than does the early- 
type population itself (Figure 7). This might indicate that 
the formation histories of BCGs involve less recent star for- 
mation than occurs in early-type galaxies of similar mass; 
below average recent star formation rates would be con- 
sistent with dry merger formation histories. 

BCGs define the same g — r color-magnitude relation as 
the bulk of the early- type population (Figure 8). Hence, if 
BCGs form from dry mergers, then BCG progenitors must 
lie redward of this relation (else the product of the merger 
would lie blueward of this relation) . In turn, this suggests 
that BCG progenitors host older stellar populations than 
is typical for their luminosities. This may provide an inter- 
esting constraint on dry merger models. To illustrate why, 
the Appendix compares our measurements with the color- 
magnitude relation predicted by semi-analytic galaxy for- 
mation models. The models suggest that BCGs should lie 
slightly blueward of the relation defined by the bulk of the 
population (Figure 12) — an effect which is not seen in the 
data (Figure 8). 

The C4 BCGs also define a tight correlation between lu- 
minosity and dynamical mass oc R e a 2 , although this cor- 
relation is slightly steeper than it is for the bulk of the 
early- type population (Figure 5). This, with the steeper 
R e — L relation, implies that they define a shallower a — L 
relation than the bulk of the population (Figure 6). These 
scaling relations are altered even for objects which are well 
described by classical deVaucouleur profiles (Figures 4- 
10). These changes to the scaling relations are consis- 
tent with previous work based on smaller samples of local 
BCGs (Malumuth & Kirshner 1981, 1985; Oegerle & Hoes- 
sel 1991). 

The change in slope of the a — L relation may be related 
to studies of curvature in the M, — a relation. This is be- 
cause the most massive galaxies are expected to host the 
most massive black holes. This expectation is based on 
the strong and tight correlation between M, and the ve- 
locity dispersion a of the spheroid/bulge which surrounds 
it: roughly M, oc a 4 (e.g. Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Tremaine et al. 2002). M. also 
correlates with the luminosity of the bulge, although this 
correlation is not as tight (Magorrian et al. 1995). The 
curvature we see in the a — L relation (Figure 6), suggests 
that a single power law may be a better description of the 
M, — ibuigo relation than of M, — a. So it is interesting 
that Ferrarese et al. (2006a, b) find that the mass of the 
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central compact object correlates better with bulge mass 
than with a or L. 
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APPENDIX 

A: COMPARISON WITH SEMI-ANALYTIC GALAXY 
FORMATION MODELS 

The introduction states that semi-analytic galaxy for- 
mation models make strong assumptions about the for- 
mation histories of BCGs. This Appendix shows that 
although the models produce massive galaxies which lie 
on a well-defined color magnitude relation, different mod- 
els disagree on the slope and zero-point of this relation 
(Figure 11). Nevertheless, the models do agree that the 
very most massive objects lie slightly blueward of the 
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Fig. 11. — Color-magnitude relations defined by massive 
galaxies in the Durham and Munich semi-analytic galaxy 
formation models. Crosses show the range spanned by 
the second and third quartiles in color as a function of 
magnitude (the horizontal error bar shows the width of 
the magnitude bin). 
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Fig. 12. — Predicted g — r color-magnitude relation in 
the Munich model. Crosses show the range spanned by 
the second and third quartiles in color as a function of 
magnitude (the horizontal error bar shows the width of 
the magnitude bin). Solid line shows the relation in the 
SDSS. The average color is redder than in the SDSS, and 
the trend for the most massive objects to be slightly bluer 
is not seen in the SDSS. 



color-magnitude relation defined by the bulk of the bulge- 
dominated population. The SDSS early-type BCGs stud- 
ied in the main text do not show such a blueward trend 
(Figure 8). 

The 'Munich' and 'Durham' models we study here were 
kindly made available by Croton et al. (2006) and Bower 
et al. (2006) respectively. Both are based on the same un- 
derlying dark matter distribution — that of the Millennium 
Simulation (Springcl et al. 2005). For the Durham models, 
we selected BCGs of halos more massive than 10 14 /i _1 M Q ; 
this resulted in about 3200 objects from a (500/i _1 Mpc) 3 
volume. The color magnitude relation defined by these 
objects is shown in the left-hand panels of Figure 11. Se- 
lecting BCGs from the Munich models is less straightfor- 
ward, because information about halo masses is not pro- 
vided. However, the stellar masses of the Durham model 
BCGs are all greater than 1.6 x 10 10 /i _1 Mq, and so the 
right hand panels of Figure 11 show the color-magnitude 
relation of the objects which, in the Munich models, have 
stellar masses greater than 2 x 10 10 ft. -1 MQ. The difference 
between the Durham and Munich models is striking, and is 
not particularly sensitive to this cut. This illustrates that 
the semi-analytic modelers have not yet converged on an 
unambiguous prediction for the slope and zero-point of the 
color-magnitude relation. 

For a more direct comparison with the results presented 
in the main text, Figure 12 shows the predicted g — r color 
as a function of r-band magnitude in the Munich models. 
Comparison with Figure 8 shows that the predicted colors 



are redder than in the SDSS. In addition, the trend for 
the most massive objects to be slightly bluer is not seen in 
the SDSS. If color gradients are responsible for the overall 
offset in color, then they must change along the color- 
magnitude relation to reconcile the predicted bluer colors 
of the most massive objects with the observations. 

B: BCGS AND THE MOST LUMINOUS GALAXIES 

So far we have focussed on the fact that BCGs appear 
to have larger than expected sizes for their luminosities. 
It is, of course, possible that a single power-law is not a 
good description of the R e — L relation of normal early- 
types. Figure B13 presents a study of this correlation for 
early- types which are not BCGs. The top left panel shows 
the joint distribution of R e and L in the full Bernardi et 
al. (2003a) sample, using the photometric reductions of 
Hyde et al. (2006). Notice that the lower R e envelope 
of the distribution appears to curve upward at large i, 
whereas the upper R e envelope is approximately constant. 
The solid line shows the result of fitting a single power to 
(R e \L). Upper (magenta) and lower (cyan) jagged lines 
show the median size in small bins in L in the C4 BCG 
sample and in the total sample. At large L, both samples 
curve upwards, away from the single power law fit. 

The upwards curvature in the full sample (normal plus 
BCGs) is qualitatively consistent with the curvature in 
the isophotal size-luminosity relation for cluster ellipticals 
recently noted by Cypriano et al. (2006). (In principle, 
curvature in the isophotal size-luminosity relation can re- 
sult even if there is no curvature in the R e — L relation; we 
have checked that this cannot account for the curvature 
they report.) What causes this curvature? 

The R e — L relation is a number-weighted average of 
the normal and BCG relations. So, if the fraction of ob- 
jects which are BCGs increases with increasing luminos- 
ity, then this will produce curvature in the same sense 
as seen in Figure B13. If all of the curvature is to 
come from this effect, then the fraction of BCGs must 
increase with luminosity approximately as /bcg(-^i-) = 
[1 + exp(M r + 23.5)] -1 ; about half of all objects with 
M r = -23.5 must be BCGs. 

To illustrate that a substantial fraction of the most lumi- 
nous galaxies are BCGs, Figure B14 shows 1.5' x 1.5' fields 
centered on a random selection of objects with M r < —23 
in the SDSS/Bcrnardi et al. sample. The fields tend to 
be very crowded, with many objects in the field having 
the same color as the more luminous central object, thus 
strongly suggesting that the central object is a BCG. 

Unfortunately, the vast majority of objects with such 
luminosities in the SDSS lie at redshifts which are sub- 
stantially larger than those over which the C4 catalog is 
complete — all the objects in Figure B14 are at z > 0.2; 
the C4 catalog is only complete to about z ~ 0.12. This 
makes it difficult to quantify how the BCG fraction in- 
creases with luminosity. Nevertheless, the top right and 
bottom left panels in Figure B13 show an attempt to ad- 
dress this issue. The points in the top right panel show the 
result of limiting the Bernardi et al. sample to the same 
redshift range as the C4 catalog (z < 0.12), and then re- 
moving the C4 BCGs as well as the objects identified as 
BCGs in the maxBCG catalog of Koester et al. (2007). 
The solid line is a single power-law fit to the points in this 
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panel; it is essentially the same as in the previous panel. 
The jagged magenta line shows the same BCG relation as 
in the previous panel. Although the lower envelope still 
curves upwards, the upper envelope is no longer as con- 
stant as before — removing the BCGs has depleted this up- 
per envelope substantially. Nevertheless, the cyan jagged 
line suggests that the objects which survived the BCG cuts 
do still lie above the power-law relation at M r < —23. 
What are these non-BCG objects with abnormally large 
sizes? 

Recall that the C4 catalog only contains clusters with 
more than ten luminous galaxies, so the early-type sample 
which remains still contains BCGs of less massive clus- 
ters and groups. Since there is no reason to believe that 
the formation histories of the BCGs of massive clusters 
are particularly different from those of rich groups, it is 
possible that such objects contribute to the curvature in 
the top right panel of Figure B13. However, a visual in- 
spection of the images of the objects with M r < —22.5 
showed that many are actually rather complex. Although 
some have multiple near-neighbours of similar color (pre- 
sumably these are the lower mass groups which are not 
included in the C4 or maxBCG catalogs), others appear 
to have double nuclei, or other unusual features. Hyde 
et al.'s estimates of the sizes and luminosities do not ac- 
count for such complexities, making their estimates for 
these quantities suspect. Therefore, on the basis of these 
images, we classified objects as having suspect photome- 
try or not. Figure B15 shows 1.5' x 1.5' fields centered on 
the objects with M r < —23 which we removed, and Fig- 
ure B16 shows similar fields centered on the objects which 
we kept. 

The jagged cyan line in the bottom left panel of Fig- 
ure B13 shows the R e —L relation of the objects which we 
felt had more reliable photometric reductions. Notice that 
the difference between the jagged and solid black lines, 
which was obvious in the top right panel, has gone away, 
even though the actual values of R e and L played no role 
in determining which objects were kept. 

As a check, the bottom right panel of Figure B13 shows 
the SDSS reductions of the same objects shown in the 
top right panel. Recall that these reductions slightly un- 
derestimate the luminosities and sizes of bright objects in 
crowded fields (such as BCGs); the solid line which shows 
a fit to these points has a shallower slope than in the other 
panels. The BCGs clearly lie above this relation at large 
L, whereas the non-BCGs do not. 

Thus, although our results do not exclude the possiblity 
that there is real curvature in the R e —L relation of normal 
early-type galaxies, they do suggest that to estimate this 
curvature reliably, one must be careful to ensure that the 
sample contains no BCGs or contaminated galaxies. In 
principle, one can use recent halo model interpretations 
(see Cooray & Sheth 2002 for a review) of the luminosity 
dependence of galaxy clustering to estimate this effect, but 
this is beyond the scope of the present work. 
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Fig. B13. — Comparison of the R e — L relation denned by BCGs and the rest of the early type galaxy population. Magenta 
jagged solid line (same in all but bottom right panel) shows the median BCG size as a function of luminosity for bins in 
L r of width 0.2 mags. Jagged cyan line in each panel shows a similar analysis of the full sample of early types (top left; 
BCG + non-BCG), of a subsample limited to z < 0.12 in which C4 BCGs plus max-BCGs have been removed (top right) 
and of the same subsample when objects with suspect photometry have also been removed (bottom left). Straight solid 
line shows a power-law fit to the points in each panel: it is essentially unchanged (and well described by equation 1) in all 
but the bottom right panel. Bottom right panel is similar to top right, but with SDSS parameters instead of those from 
Hyde et al. (2006). 



14 



0. 24827-23. 83/1. 22 
J17284S. 11 + 610730. 2 


0.2821/-23.90/1.27 
J112916. 81-004244.4 


0.2041/-23.95/1.17 
J094434. 45-0 04429. 9 


0.2118/-23.83/1.17 
J081045. 58+493540. 9 


0.2514/-23.78/1.14 
J171616. 79 + 590731 


• • 
• 

• 


• 


• 


V 


1 


0. 20 35/-Z4. 01/1.32 
J110806. 11 + 001653. 5 


0. 2493/ 23. 92/1. 36 
J00423S. 09 + 153104. 5 


0.2238/-23. 74/1.16 
J074255. 64 + 390310 


0,1680/-23.80/1.19 
J125754. 3-015851. 7 


0.1772/-23.B1/1.21 
J082632. 84+495738.4 


V 

% 

■ 


% 

■ 


• 




( 


0.2621/-23.81/1.10 
J141147. 04 + 002420. 2 


0.2001/-23.77/1.16 
J075113. 87 + 383548. 3 


0.2459/-23.95/1.36 
J015248. 27 + 150228. 2 


0.2717/-23.78/1.11 
J095050. 05-005315. 8 


0.1800/-24.01/1.20 
J101534. 39-005251. 9 


• 

■ 

• 

• 


• 




• 

• * • 


» 



Fig. B14. — 1.5' x 1.5' fields centered on a random selection of objects with M r < —23 in the SDSS/Bernardi et al. 
(2003a) sample. The fields tend to be very crowded, with many objects in the field having the same color as the brighter 
central object, as expected if the central object is a BCG. Band at top of each image shows the redshift, the absolute 
magnitude in r-band, log 10 R e of the object and its SDSS ID. 
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Fig. B15. — Objects with M r < —23 and z < 0.12 in the SDSS which are in crowded fields or are contaminated/peculiar 
galaxies: these objects were present in the non-BCG sample shown in the top right panel of Figure B13, but were removed 
to make the bottom left panel of that Figure. Images have been ordered by luminosity: the band at the top of each image 
gives the absolute magnitude in r-band, log 10 -^e an d the object ID. 
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Fig. B16. — Objects with M r < —23 and z < 0.12 in the SDSS which are in relatively clean fields or do show strong 
contamination: these objects make the non-BCG sample in both the top right and bottom left panels of Figure B13. 
Images have been ordered by luminosity: the band at the top of each image gives the absolute magnitude in r-band, 
log 10 R e and the object ID. 



